Algebro-geometric approach in the theory of integrable 

hydrodynamic type systems 



Maxim V. Pavlov 

Abstract 

The algebro-geometric approach for integrabihty of semi-Hamiltonian hydrody- 
namic type systems is presented. This method is significantly simplified for so-called 
symmetric hydrodynamic type systems. Plenty interesting and physically motivated 
examples are investigated. 
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1 Introduction 

The theory of the integrable hydrodynamic type systems 



was estabhshed by S.P. Novikov, B.A. Dubrovin (see [6]) and S.P. Tsarev (see [31]). 
This differential-geometric approach has been developed also by E.V. Ferapontov (see, 
for instance, [7], [10] and other references therein), I.M. Krichever (see, for instance, [18] 
and other references therein), O.I. Mokhov (see, for instance, [23] and other references 
therein) and by the author (see, for instance, [28] and other references therein). Also, B.A. 
Dubrovin and I.M. Krichever (see, for instance, [5], [19] and other references therein), 
Yu. Kodama and J. Gibbons (see, for instance, [14], [16] and other references therein) 
used algebro-geometric approach in the theory of integrable hydrodynamic type systems, 
which are dispersionless limits of integrable dispersive equations (or in most general case 
the Whitham equations, i.e. obtained by the Whitham averaging method of multi-phase 
solutions of dispersive systems). Thus, all corresponding information like the Riemann 
surfaces, a quasi-momentum and a quasi-energy can be reconstructed. Then (as usual in 
the algebro-geometric approach) generating functions of conservation laws and commuting 
flows can be found automatically. 

This paper is devoted to the algebro-geometric approach for hydrodynamic type sys- 
tems, whose origin is unknown. For simplicity in this paper we restrict our consideration 
on symmetric hydrodynamic type systems, because just in this case a generating func- 
tion of conservation laws in fact is given in advance. In all other cases derivation of 
a generating function of conservation laws is separate and complicated computational 
problem, which will be investigated in details in another publication devoted integrable 
hydrodynamic chains. The next step is a computation of the equation of the Riemann 
surface. The corresponding linear ODE system can be solved if it is invariant under a 
some Lie group symmetry. For instance, if the symmetric hydrodynamic type system is 
homogeneous (this is a typical formulation of physically motivated examples), then the 
corresponding generating function of conservation laws and the equation of the Riemann 
surface are homogeneous too. 

Moreover, the generalized hodograph method established by S.P. Tsarev in [31] is 
based on a concept of the Riemann invariants. In this paper we suggest an alternative 
approach based on a conservative form of hydrodynamic type systems. Most physically 
motivated problems are given in such form. 

The paper is organized in the following order. In the second section a semi-Hamiltonian 
(integrability) property for hydrodynamic type systems is reformulated for a conservative 
form. The method allowing immediately to construct generating function of conserva- 
tion laws and the corresponding Riemann surface is established. In the third section we 





(1) 
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briefly describe Tsarev's observations improving plenty calculations. In the fourth sec- 
tion the chromatography system as the most interesting and very complicated example 
of symmetric hydrodynamic type systems is investigated. In the fifth section the general- 
ized hodograph method adopted to a conservative form is considered in details. Several 
different sub-classes of hydrodynamic type systems are presented. In the sixth section 
homogeneous hydrodynamic type systems are considered. In such case a computation 
of the Riemann surface can be found in quadratures. In the seventh section integrable 
hydrodynamic chains as a natural generalization of the symmetric hydrodynamic type sys- 
tems are discussed. In the eighth section the Hamiltonian chromatography hydrodynamic 
type system is considered. Corresponding hydrodynamic chain is found as well its local 
Hamiltonian structure. In the ninth section the linearly degenerate case is investigated. 
In the conclusion the main open problem is considered. 



2 Symmetric hydrodynamic type systems 

The hydrodynamic type system (1) written in the Riemann invariants 

ri = /iXrK, 2 = l,2,...,iV (2) 



d.-T^ = d.-i^.^ ^^J^k (3) 



is integrable by the generalized hodograph method (see [31]) iff the semi-Hamiltonian 
condition 

-r : = Ok — ■ 

fi" — yU* fL^ — fi' 

is valid identically (here dk = d/dr^\ the semi-Hamiltonian condition also can be written 
in arbitrary field variables u^, see [27]). Then the hydrodynamic type system (1) has 
infinitely many conservation laws and commuting flows 

rl = w\v)rl, (4) 

parameterized by arbitrary functions of a single variable. 

Remark: Characteristic velocities w*(r) satisfy the linear PDE system (see [31]) 

dkw^ = ^^^{w^-w^), i^K (5) 

which cannot be solved explicitly in general case, because this is a linear system with 
variable coefficients. However, in this paper we are able to avoid this problem, because 
the generalized hodograph method can be formulated via arbitrary field variables (see 
the beginning of the section 6). The algebro-geometric approach is based on the concept 
of the Riemann surface, where all conservation laws and commuting flows can be found 
(see below) via the Riemann invariants as well as conservation law densities u^, which 
appear in this framework in a natural way. The system (5) has the general solution 
parameterized by N arbitrary functions of a single variable. The generalized hodograph 
method established by S.P. Tsarev (see [31]) leads to the general solution written in 
implicit form as an algebraic system 

x + nH = w\ i = l,2,...,Ar (6) 
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for the semi-Hamiltonian hydrodynamic type system (2). In this paper several tools useful 
for constructing of commuting flows are suggested in the section 5. 

Plenty hydrodynamic type systems are known from physical applications, which have 
the very special conservative form 

Mi = a,V^(n^n^...,n^;nO, (7) 

where the sole function ip^u^^u^, ...,u^;p) is invariant under permutation of the first 
entries (also we restrict our consideration on the case, when the function ?/'(u;p) is 
nonlinear with respect to p, see details in the last section 9) . We call such hydrodynamic 
type systems the symmetric hydrodynamic type systems. The theory presented below can 
be easily extended on more general symmetric classes, for instance on 



V 



where N and M are arbitrary integers. Moreover, we shall demonstrate below (in a set 
of examples) that this is not necessary restriction. However, in the symmetric case this 
theory is very effective. 

Let us introduce the matrix given by 

dip dil)\ dijj 



and formulate the phenomenological algebro-geometric approach. 

Statement 1: If the symmetric hydrodynamic type system (7) is integrable, then this 
system has the generating function of conservation laws 

Pt = d^ip{u^,u'^,...,u'^;p). (9) 

If the generating function of conservation laws (9) is consistent with the hydrodynamic 
type system (7), then 

where the matrix B^{u;p) is an inverse matrix to the matrix Af (see (8)). This is 
ODE's of the first order for every fixed index i, where any of them can be written in the 
form dy/dx = f{x,y). 

Let us introduce the function X{u^, u"^, u^\p) determined by N linear PDE's of the 
first order (see (10)) 

d\ ^ dip dX _ ^ 
* du'' du^ dp 

Definition 1: The function A(u;p) satisfying (11) is said to he the equation of the 
Riemann surface. 

If the linear PDE system (11) has an integration factor then the equation of the 
Riemann surface can be found in quadratures 



dX = — 
op 



dp-BJ^-^du'' 



4 



If, for instance, the hydrodynamic type system (7) is homogeneous (see examples below), 
then functions ■?/'(u;p) and A(u;p) are homogeneous too. Then the integration factor can 
be found by using the Euler theorem 

dX ^dX 



X = Ptt- + w 



dp du'^ 

Statement 2: A deformation of the Riemann surface determined by the equation 
X{u]p) satisfies the Gibbons equation 

X,-^X, = ^[p,-dM^-p)]. (12) 

The Gibbons equation (first introduce in [11]) has three distinguish features: 

1. if one fixes A = const (free parameter), then one obtains (9), 

2. if one fixes p = const (free parameter), then one obtains the kinetic equation (a 
coUisionless Vlasov equation) written in so-called Lax form 

X - a H} - 

' dx dp dp dx ' 

where H = iplu^p). 

3. if one choose coordinates, which are the Riemann invariants r* [i = 1,2,.. .,N) 
determined by the condition dX/dp = (see (11)), then the corresponding hydrodynamic 
type system (7) can be written in the diagonal form (2) 

ri = ^\p=p^rl, 2 = 1,2,. ..,7V, (13) 

where the corresponding values can be expressed via these Riemann invariants r^. In 
this algebro-geometric construction the Riemann invariants are the branch points r* = 
Md\/dp=o of the Riemann surface (exactly as it is in the Whitham theory, see [5] and 
[19]). 

Remark: The characteristic velocities //^ of hydrodynamic type system (1) can be 
found from algebraic system 

det\vUu)-p6l\=0. (14) 

All of them must be distinct in agreement with Tsarev's assumptions (see [31]). Thus, 
if the hydrodynamic type system (7) is integrable by the generalized hodograph method, 
then the values are given by (see (13)) 

di/j 

where characteristic velocities are determined from (14) 

detAf(u;p) = 0. (16) 
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Suppose the hydrodynamic type system (7) is semi-Hamiltonian. Then such system 
must have N series of conservation laws, which can be obtained by the substitution of the 
formal series 



in (9). The compatibility conditions of the first TV extra conservation laws 



(17) 



dtvUn) = 



V u 



dip , 
dp 



\p=u' 



with the hydrodynamic type system (7) are equivalent the semi-Hamiltonian property. 

Main statement: The symmetric hydrodynamic type system (7) is semi-Hamiltonian 
iff the compatibility condition di{dkp) = dk{dip) is fulfilled. 

Comment: Computation of this compatibility condition can be made in the coordi- 
nates (see (7)) or in the Riemann invariants (see (13) and details in the Conclusion). 
In both cases the nonlinear PDE system in involution coincides with integrability crite- 
rion for hydrodynamic type systems following from existence of conservation laws and 
vanishing Haantjes tensor (see [15], [27]). 



3 Tsarev's observations 

Let us consider the dispersionless limit of the vector NLS (see [34]) 



i\2 



7f^ = d^iuY), 



l,2,...,iV. 



(1^ 



The eigenvalue-eigenfunction problem (cf. (16)) is 



(m* - n)6ik 1 



Thus, 



— 



/ ^ . (19) 

n — w — ' (/i — uf-y ^ 

The first Tsarev observation [32] is that the sum of the last equations yields an 
expression determining characteristic velocities fi'' (see (14) and (16)) via very compact 
formula 

i = y . 

^ (/i - u'^y 

The second Tsarev observation [32] is that the above expression can be integrated 
once with respect to 



(20) 

Comment: This equation of the Riemann surface can be obtained directly from 
the spectral problem for the vector NLS (see [34]). The characteristic velocities fi^ are 
determined by the condition dX/dfi = 0, where the Riemann invariants r* = \\d\/d^=o = 
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X\^=^i are branch points of the Riemann surface (see [11])- A deformation of the Riemann 
surface is described by the Gibbons equation (see [11] again) 



(21) 



which connects (18) with (20). 

The third Tsarev observation [32] is that the fiat diagonal metric of the hydrodynamic 
type system (18) written in the Riemann invariants (2) 



is given by (see also [31]) 



i = 1,2, ...,2N 



gu = res — d^i. 

\=ri \Ojl J 



(22) 



Indeed, since the Hamiltonian structure of the dispersionless limit of the vector NLS is 

i 9h .ah 



the diagonal metrics (in Riemann invariants) 



is given by 

where (see (20) and cf. (19)) 

(9r* 



N k 

9" = 2E^ " 



Thus, = 9^A/9/i^|^=^i in agreement with (22). 



4 Generalized chromatography 



The chromatography process (see, for instance, [10]) is described by the hydrodynamic 
type system 



2 = 1,2, 



(23) 



where a,(3,e and 7^ are constants. All results presented in this section generalize results 
from [10] (see formulas 1, 3, 4, 5). This hydrodynamic type system has the obvious 
couple of conservation laws 



dt 



fc\/3-o+l 



(/?-« + 1)9, 

7 



/3{a + e - pe) i^P^ ^ 



U 



n\a+f3—l 



a + p-l 

where A = 1 + ^ '^f.{u^Y . 

If a = /3 — 1, then the hydrodynamic type system (23) is HamUtonian (see [6], [28]) 



-d — 

7j ^ du" ' 



(24) 



where the momentum density is Y^^^Jyvf^Y' ^^"^ ^^e Hamiltonian density is h = A^~^. 
Suppose the above hydrodynamic type system is integrable for some values of constants 
e and 7^. Then the generating function of conservation laws given by 



Vt = d, 



should exist. It is easy to check that the compatibility conditions di{dkp) = dk{dip) are 
valid iff a = Pe, where (see (10)) 



dp 



7i(«*)^~V 



n\a+l3-l 



(25) 



All first derivatives (see (11)) 
— = ^{n,p)p^^ 



dX 

^ = V^(u,p) 



n\l3 



'y^n'^/Se—l pPe—l 



are determined up to integration factor (y9(u,p), which is not found yet. Then the equation 
of the Riemann surface A(u,p) can be found in quadratures 



d\ = Lp{\i,p) dp 



,fl+i 



(3e 



{w'^)w=i-^dw' 
u;" - 1 



if the integration factor is ^{p) = p ^ ^ (here we use the substitutions p = qi^^-^ and 



= {w'^q)!^^-^). Then, the integrable hydrodynamic type system (23) 

,i\f3e 



i _ \U, 

is connected with the equation of the Riemann surface 

-/3 1 



1,2, ...,iV 



(26) 



A 



P 



/5 p 



where 2Fi{a, b, c, z) is a hyper-geometric function and S = P/{Ps — 1). If 5 = m/n, where 
m and n are integers, then the equation of the Riemann surface A(u, p) can be found in 
elementary functions. 
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Remark: If /9 — oo, then the hydrodynamic type system (26) reduces to 

su'' 

The corresponding equation of the Riemann surface is 

where p = Ing and m" = ln(w"g). Moreover, if £ = 1, the above hydrodynamic system 
has the local Hamiltonian structure (24) with the Hamiltonian density h = InA. The 
above equation of the Riemann surface reduces to 

A = e-^-5^7„ln(e""-^-l). (28) 

Remark: If £ ^ 0, then the hydrodynamic type system (26) reduces to 

«I = 5.1n^^^|^, ^ = l,2,...,iV, 

where the constants 7^ are removed by appropriate scaling of the field variables m*^. The 
corresponding equation of the Riemann surface is 

^ = ^ - $^(^")'F(1, -A 1 - /3, 

where = w'^p. If f3 = 1, the corresponding hydrodynamic type system (26) has the 
local Hamiltonian structure (24) with the Hamiltonian density 

h = ^M"*(lnM'" - 1) - (1 + ^M'")[ln(l + XI"") - 



5 The generalized hodograph method 

If the hydrodynamic type system (1) is semi-Hamiltonian, then the general solution pa- 
rameterized by arbitrary functions of a single variable is given in an implicit form by 
the algebraic system (cf. (6); see [31]) 

x5l + tvl{u) = wi{u), (29) 

where wl{u) are characteristic velocities of an arbitrary commuting flow. 

However, in this paper above, we present the approach producing the generating func- 
tion of conservation laws only. Thus, we need to extend this mechanism to produce the 
generating function of conservation laws and commuting flows simultaneously. In this pa- 
per we restrict our consideration on three sub-cases connected with the Egorov conjugate 
curvilinear coordinate nets (see [28]), with the orthogonal coordinate nets (see [31]) and 
with so-called "mirrored" conjugate nets (see below). The general case will be considered 
elsewhere. 
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Nevertheless, the first step is a description of N series of conservation laws (see (17)). 
They can be obtained by expansion in the Biirmann-Lagrange series (see, for instance, 
[22]) at the vicinity of each singular point. 

Theorem 1 [22]: The analytic function 

y = yi{x- Xq) + y2{x - x^f + ?/3(x - x^f + ... 
can he inverted {y{x) — * x{y)) as the Biirmann-Lagrange series 



X = xo + xiy + x^y^ + x^^y^ + 



whose coefficients are 



1 dJ'-~'^ f X - xo\"' 
^" = rTiim3— T— , n = l,2,... (30) 



nlx^xodx"^"^ 

For example, the so-called "waterbag" hydrodynamic type system (see [14], [16]) 



< = d.i^^ + J2'ka'j (31) 

is connected with the equation of the Riemann surface 

>^ = p-J2ekHp-a''). (32) 

The main (so-called "Kruskal") series of conservation law densities can be obtain by 
substitution of the Taylor series 

Hg Hi H2 

p = X---^-^-... (33) 

into the above expression if Se^ = 0. If Se^ 7^ 0, then at first, the above equation of the 
Riemann surface must be replaced on 

A-^£felnA =p- J]efcln(p-a^), (34) 

because the Gibbons equation is invariant under scaling A A (A). In both cases are 
polynomials with respect to field variables a'^. 

Also, the "waterbag" hydrodynamic type system has N infinite series of conservation 
laws. At first let us rewrite the equation of the Riemann surface in the form 



A = {p - a')e-P/''Yl{p 



for any fixed index i. Then the infinite series of conservation laws (17) 

p« = a' + hf (a) A + hf (a) A' + 4' V) A' + . . . (35) 
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can be obtained with the aid of Biirmann-Lagrange series (see [22]), which coefficients 
are determined by (see (30)) 



1 



Thus, the first conservation laws are 



, . . . 



„2a7ei / \ 

/i« = e«V.. _ a'^)-'^/'', hf = 1 - 5^-^ n^^' - a'=)-2-^/^' 

5.1 Quasi-symmetric form 

The dispersionless hmit of the vector NLS (18) is a degeneration of the "waterbag" 
hydrodynamic type system (31). A substitution of the expansion (see, for instance, [2]) 
= u^ + Tj^/e^ + ... in (32) 

N 

A = /i - yek In ^ 

k=i 

yields (20) if oo. 

Let us consider again the generating function of conservation laws (see (21)) 

/i, = 5x + A'^ . (36) 

for the Benney hydrodynamic chain (see [1]) 

A^t = A^/^ + kA^-'Al fc = 0,1,2,... (37) 

and substitute Taylor series (35) 

^(^ = a' + Xb' + X^d + ... (38) 

1. Suppose the function depends on field variables only. Then corresponding 
hydrodynamic type system is 

However, this hydrodynamic type system is integrable iff the function A^{a.) satisfies some 
nonlinear PDE system 



fa* - a'' 
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which is a consequence of the compatibihty conditions di{dkfi) = dk{difi), where 

jjL — /i — a" ^ 

Several such particular choices are described above. 

2. Suppose the function depends on field variables a*^ and M field variables 
6^ (where M must be not exceed A^) only. Then corresponding hydrodynamic type 
system is 



a. (^^ + A°(a,b)^, z = l,2,...,iV, 



^ = d^ia^V), j = l,2,...,M. 

Simplest such example is (18). This procedure can be extended on other auxiliary field 
variables from (38). For instance, the third such sub-case is 

4 = 5. (^ + A°(a,b,c)^, ^ = 1,2,..., AT, 
bi = 9,(aV), j = l,2,...,M, 



aV + -(6' 



k\2 



fc = l,2,...,i^. 



where K ^N. 

3. Suppose the function A^ (see (37)) is a function of the field variables only. Then 
corresponding hydrodynamic type system is 



A', = d,A\a.) 



9. 



A' 



1,2,...,N. 



(39) 



Suppose the function A is a function of the field variables a only. Then corresponding 



hydrodynamic type system is 



A° = dxA\ 



Al = AHsi) + 



{A'f 



al = 



J\2 



A' 



1,2,..., AT. 



Thus, in general case one can consider M first moments A^ and K first sets a*^, 6*^, c*^, 
... (where the index run all values from 1 up to A^^, A; = 1, 2, ...,K) as field variables 
of corresponding hydrodynamic type systems. In such case just the latest moment A^ 
depends on sets a*^, 6*^, c*^, ... 

Nevertheless, the theory established above is still working for these hydrodynamic type 
systems, because all of them have the same generating function of conservation laws (36). 

Example: The dispersionless limit of the vector Yajima-Oikawa system (see [28]) is 
the hydrodynamic type system (cf. (18) and (39)) 

,k\2 



A' 



ril = dxiu^v ) 



k = 1,2,..., AT. 



(40) 
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The eigenvalue-eigenfunction problem is 



-/i 1 

1 (n' - fi)6ik 

ri'6ik (u* - fi)Sik 



r 



Thus, 



r] 



where 2N + 1 eigenvalues are solutions of the discriminant equation 



In accordance with the Tsarev observations, let us integrate this expression once. Indeed, 



7]" 



/i — 



is the equation of the Riemann surface connected with the hydrodynamic type system 
(40). The factorized form of this equation 



N+l 

n (/i - a-) 

k \ n=l 

N 

n (/^ - ^^') 
k=l 



(41) 



yields 27V + 1 independent series of conservation laws 



(k) 



u'^ + A?7''(u, a) + \w^{\i, a) + 



fc = l,2,...,iV, 



in) 



a" + A6"(u, a) + Ac"(u, a) + n = 1, 2, iV + 1, 



whose coefficients can be found with the aid of the Biirmann-Lagrange expansion (see 
(30)). 

Remark: While the conservation law densities r]'^ are coefficients of the first series, 
the function is a coefficient of the Kruskal series for (41) at infinity A — > oo, /i ^ cxd 



A = ^ + A° + - 
2 /i 



+ ..., 



where = Sr/" (see (39) and (40)). 



5.2 Hamiltonian formalism 

The Hamiltonian formalism of hydrodynamic type systems was established in [G] for a 
local case, and developed later for a nonlocal case in [7], [9]. The generating function of 
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conservation laws for the symmetric hydrodynamic type systems (7) is given by (9). Sup- 
pose for simplicity without lost of generality that the field variables are flat coordinates. 
Then (7) can be written in the symmetric Hamiltonian form 



9h 



where the Hamiltonian density h is a symmetric expression under an index permutation 
{si and 7j are some constants). Thus, the generating function of commuting flows is given 
by 

i ( dp sr^ dp 

Then infinitely many particular solutions are given by the generalized hodograph method 
(29) 



j d I dh dh 



oo N 

EE 

71 = 1 S = l 



d / dh); 



is) 



du'^ \ du^ 



dh 



du^- 



where agn are arbitrary constants. 

Example: The Hamiltonian exponential chromatography (27) 



ul = d^ 



i = 1,2, ...,N 



■1 + E7fce"'^ 

is connected with the equation of the Riemann surface (28) by the Gibbons equation 



A ■ 



dX 

dp 



Pt - d^ 



A 



The Kruskal series of conservation law densities can be found by the application of the 
Biirmann-Lagrange series (30) at the vicinity X —>■ 0,q —>■ 0, where q = exp{—p). In this 
case the coefficients g„ of the inverse series 



are determined by 



1 , 

= -rlim- 



n = 1,2,... 



nlq^odq"- ^ \g — ^ 7^. ln(l — ge" 
Then the Kruskal conservation law densities can be found from p = — In g 

p = p + In X = - ln[l + gsA + qsX^ + QaX"^ + ...]. 

For instance, pi = In A, p2 = A~^E7„e^"". 
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The Kruskal series can be extended in the opposite direction A — >■ oo, q oo. At first, 
the equation of the Riemann surface (28) must be replaced by (cf. (34)) 



where s = exp(— p) — S7^u'". Then the Kruskal conservation law densities p^k can be 
found exactly as for the Benney hydrodynamic chain (37) (see (34)) by substitution (cf. 
(33)) 

, Ci C2 C3 

in the above formula. Then p_k are coefficients of the series 

(1 C C C \ 

For instance, p_i = S7^m™, p_2 = S7^e~"'" + (S7^m'^)^/2. 

Any semi-Hamiltonian hydrodynamic type system has series of conservation laws. 
The equation of the Riemann surface (28) can be written in the form 



A = (e-f - e-"') exp - — + ln(e"'-P - 1) 



p ^ u . 



Then the coefficients qt of the Biirmann-Lagrange series are given by (30) 

1 rf"-^ (l-^x 1 o 

qn = —. lim . "TT^ ' n = l,2,..., 

nlq^cxp{-u')dq"-~^ \ X{q) 



where q = exp(— p). Finally, A^ series of conservation law densities p*-*-* can be found from 
the series 

p{i) ^ p{i) + inA = -ln[gi + ggA + qsX^ + q^X"^ + ...]. 
For instance, the first A^ nontrivial conservation law densities are 

p« = e-"'-5^7.1n(e"'=-'^'-l). (42) 
Commuting flows of this Hamiltonian chromatography system are 



,1 



*jfc — dx " , A; — 0, ±1, ±2, 

7j C/M* 

dhk,n 

7j du^ 

where are Kruskal conservation law densities and hk^n are A^ series of conservation law 
densities. The flrst A^ commuting flows 



ulk,n = — (9^.— -y^, A; = l,2,..., n = 1,2, A^, 



are given by (42) 



Ik \, -e" / e'^ - 
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5.3 Mirrored curvilinear conjugate coordinate systems 

If the semi-Hamiltonian property (3) is valid, then the diagonal metric ga can be intro- 
duced by 

Jh!^ = dk\nH,, i^k. (43) 

Following G. Darboux (see [4]), let us introduce the rotation coefficients of the conjugate 
curvilinear coordinate nets 

= ^^k. (44) 

Then any solution iJ^ of the linear problem (see [31]) 

= (3,kH,, I ^ k (45) 

determines commuting flow to (2), where characteristic velocities are connected with 
Hi by the Combescure transformation = Hi/ Hi. Any solution of the conjugate linear 
problem 

di^k = Pki^i, i ^ k (46) 
determines conservation law density 

d,a = iJiHi. (47) 

Definition 2 [28]: The conjugate curvilinear coordinate net determined by symmetric 
rotation coefficients (3^^ = Pj.i is called the Egorov conjugate curvilinear coordinate net. 

Theorem 2 [28]: // the integrable hydrodynamic type system (2) has the couple of 
conservation laws 

at = bt = Or,, (48) 

then corresponding conjugate curvilinear coordinate net is Egorov. 

Remark [28]: This couple is unique for each given Lame coefficients Hk. 
Corollary [T8]: Any commuting flow (4) has similar couple of conservation laws (48) 

hy fx., 

where a is the unique potential of the Egorov metric for all commuting flows. 

Definition 3: Two conjugate curvilinear coordinate nets determined by the rotation 
coefficients and (3^^ = Pki is called mirrored conjugate curvilinear coordinate nets. 

Theorem 3: // one integrable hydrodynamic type system (2) has the conservation law 

at = bx (49) 
such that another integrable hydrodynamic type system 

rl = ~^\v)rl (50) 

has the couple of conservation laws 

ay = Cz, by = B^, (51) 
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then corresponding conjugate curvilinear coordinate nets are mirrored. 

Proof: The characteristic velocities of the hydrodynamic type system (2) are (see 
(47)) 

Hi ipiHi dia 

where H^^^ is a solution of the linear system (45). Let us consider such hydrodynamic 
type system (50), whose characteristic velocities /i* are (see (47)) 

where ipf^ is a solution of the conjugate hnear system (46). The rotation coefficients can 
be found in two steps (43) and (44). Indeed, = /^^.j. 

Remark: The construction described above is symmetric. Thus, the second conser- 
vation law of the ffist hydrodynamic type system (2) 

Ct = Cx 

is given by quadratures 

dc = J2 4^^H,dr\ dC = J2 i^[^^Hl^^dr\ 
Corollary: Any conservation law 

Py Qz 

of the hydrodynamic type system (50) determines the corresponding commuting flow in 
the conservative form 

ar = p. (52) 

or in the Riemann invariants (see(4)) 

- ' < (53) 



Hi 



(2') 

of the hydrodynamic type system (2), where H^ is some solution of the linear system 
(45), d^P = Hf~^^i and d,Q = i/f V!'^- 

Thus, if two symmetric hydrodynamic type systems are related by the above link, 
then the generating function of the second hydrodynamic type system (50) determines 
the generating function of commuting flows for the flrst hydrodynamic type system (2). 

5.4 Chromatography system 

The integrable chromatography system (26) has the couple of conservation laws 
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pe + p-1 



u" 



1. The Egorov sub-case. If (3e = —1, then this is the Egorov hydrodynamic type 
system (see (48)) 

where the potential of the Egorov metric (see [28]) is 

2. The local Hamiltonian sub-case. If /5(1 — e) = l, then chromatography system (26) 
can be written in the Hamiltonian form (24). li j3 oo, the exponential chromatography 
(27) [e = 1) also has the Hamiltonian form (24), where the Hamiltonian density is 
h = InA. 

3. The nonlocal Hamiltonian sub-case (associated with constant curvature metric, see 
[9], [24]). If /3 = 2, then the conservation law density 



is the momentum density (see [24]) of the nonlocal Hamiltonian structure 

where diagonal matrix elements gik = —li^ik (where 5ik is the Kronecker symbol) and the 
Hamiltonian density 



4. The general [mirrored) sub-case. If two hydrodynamic type systems (26) are related 
via mirrored conjugate curvilinear coordinate nets, then they must have one common 
conservation law density (see (49) and (51)). Let us prove that such conservation law 
density is 

a = Y.l,{u'Y-^^^\ (54) 
If the another chromatography system (26) 

ul = d, ^, z = l,2,...,iV (55) 

' [1 + E7.(«W 



has the same conservation law density 



then transformation u^u, j3^l3,e^e can be found from two other obvious assump- 
tions (see (49) and (51)) 

u' = iff 
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Result: 



A 1 ~ 1 3 1 



Since the second chromatography system (55) has the generating function of conser- 
vation laws 

then the generating function of commuting flows is given by (52) 
where the relationship 

P = p-^' 



(56) 



is obtained by a comparison the equations of the Riemann surface (which are equivalent) 
for both chromatography systems. 

The integrable chromatography system (26) in the Riemann invariants has the form 
(13) 

Since this hydrodynamic type system has two conservation laws 



Pt = 



at = i/3e-/3- l)d.,A- 



then one can obtain, respectively 

-^^^ a m A 



diP 



{p'f-^^dMA. 



Thus, the generating function of commuting flows (56) in the Riemann invariants (53) 



= a- 



p[(pi)Pe-l _ pf3e-l^ 

is connected with the Gibbons equation 

f3e {p')'^'-^ 



+ P-Pep{0[{p'r'-'-pf''~'iC)] 



A. 



dX 
¥(A) 



pe 



{(3e - + (3 - Pe) 



P_m 
PiO) 



F (1,^,^ + 1, 



pPe-l^Q 



where we use the notation 



a 



(3e 
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Then the generating function of commuting flows written via the physical field variables 
(cf. (26)) is 



{(3e - + P - (3e) 



u 



i\l3e-l 



Substituting the formal series dr{c_) = d^i^ +Cdt{ +C^f^4 + ••• and respectively the generating 
function of conservation law densities expanded in the series (35), one can obtain infinitely 
many generating functions of conservation laws for all commuting flows. For instance, the 
first such functions are 



Pti = dx 



F I 1,C7,(7+ 1, 



P 



Remark: Suppose that some component hydrodynamic type system contains A^ — 1 
equations (see the above generation function of conservation laws) 



t^+l ^x 

'■0 



F{l,a,a + 1, 



k\f3e-l 



Then such hydrodynamic type system is integrable if its A^th equation satisfies some 
extra conditions. This hydrodynamic type system is not symmetric like (7). However, 
the algebro-geometric approach still is valid. 

Thus, if any given hydrodynamic type system contains A^ — 1 equations 



u^t = d^iP (u; , 1, 



then one should substitute the Taylor series (35) in the generating function of conservation 
laws 

Pt = d^^ (u; ^) . 

The A^th equation will be obtained by the limit (see such example at the end of the 
sub-section "Hamiltonian formalism") 



U} = drr. 



limz/; u: 1 + 



All other computations are exactly as in the symmetric case. 



5.5 Reciprocal transformations 

The concept "reciprocal transformation" was introduced by S.A. Chaplygin (see, for in- 
stance, [29] and [30]) 

dz = A{u)dx + B{u)dt, dy = C{u.)dx + D{u)dt, 

where two arbitrary conservation laws At = and Ct = preserve the gas dynamic 
equations changing adiabatic index only. The gas dynamics is the first known example in 
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the theory of integrable hydro dynamic type systems. In this sub-section we show that the 
integrable chromatography system (26) is invariant under couple of different reciprocal 
transformations. Thus, infinite sets of such systems are related by a chain of reciprocal 
transformations, which are described below. 

1. The first such reciprocal transformation is very simple 

dz = dt, dy = dx. 

Then the integrable chromatography system (26) reduces to the chromatography system 



a 



1,2,. ..,iV, 



where 

u' = {w')^A-^/^, A = A-\ 

Thus, two chromatography systems (26) (/3,e,7^) and {1/e, 1//5, — 7^) are related by the 
transformation x ^ t. 

Remark: The first conservation law 



u 



k\l3-l3e+l 



{Pe - P - l)d.,A- 



transforms into the second conservation law 

pe 



dyA'/" = a 



A~^ 



E 



ln{W 



n\l3e+l3-l 



and vice versa. 

2. The reciprocal transformation 

dz = A^'^dx + - — — r A f 



dt. 



_l3e + l3-l 

connects the chromatography system (26) and another symmetric system 



dy = [3edt 



V 



i\l3e 



V 



V 



n\f3e+/3-l 



I3e pe + p 

where = u^A~^^^. 

However, the chromatography system (26) has the commuting flow 

,i\2-l3e 



1,2,. ..,7V, 



(57) 



+ 



u 



pe p-pe + 



{51 



It is easy to verify by the compatibility conditions dti{ul) = dt{u[i). Then the hydrody- 
namic type system (57) has the commuting flow 



dx 



1,2,. ..,iV 



[i-E7.(^T]'/^-^' 

which is the chromatography system (26) again. Thus, two chromatography systems (26) 
(/3, £, '-yj.) and {P,2/P — e, — 7^) are related by the above reciprocal transformation, where 
A = A-i = 1 - S7„(i;")^. 

Applying both reciprocal transformations iteratively, one can construct a link between 
the chromatography systems (26) with the distinct indexes P and e. 
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6 Homogeneous hydro dynamic type systems 

Another hydro dynamic type system 



arising in the chromatography (see formula 6 in [10], j3 and 7„ are arbitrary constants) is 
invariant under scahng of field variables n'^ cu^ (and appropriate scaling of independent 
variable t or x). We call such hydrodynamic type systems as homogeneous. Plenty 
physically interested hydrodynamic type systems belong to this class. 

The existence of the corresponding generating function of conservation laws 

yields (see (10)) 

dp ^ liP^ ( 1 -sp 7fc(^^T"^ / 



Also, the equation of the Riemann surface A(u,p) must be invariant under extended scaling 
cu^, p cp. Since (see (11)) 



du' (3u'[pP-^ - [u^y-^] V (3^ p^-i - {u^y-^ ) dp' 



then the equation of the Riemann surface A(u,p) can be found in quadratures 

In A = Inp + — — > In -— , 

(/? - l)(/5 + E7m) ^ pf^^-iu^^y^ 

where we used the Euler theorem 

A = pXp + '^M^Afc. 

Since the Gibbons equation (12) is invariant under point transformation A A(A), then 
the above equation of the Riemann surface A(u,p) can be written in the form 

A=/-s-n(i-^)"". 

where q = p^"^, = {u^y~^. 

6.1 The Kodama hydrodynamic type system 

The Kodama hydrodynamic type system (see [16]) 

k 



at = 



m=l 



A: = 1,2,..., AT, 
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where S^^ is the Kronecker symbol, is a homogeneous (The Euler operator is E = T.(N + 
k — 2)a^dk) hydrodynamic reduction of the Benney hydrodynamic chain (37). Thus, the 
corresponding generating function of conservation laws is (36) (the generating function of 
conservation laws for any integrable hydrodynamic chain is unique for all its hydrodynamic 
reductions). The Kodama hydrodynamic type system can be obtained from the above 
generating function by substitution of the Taylor series 



V = 

Then 



ai + Aa2 + AV + ... + A^^ia^ 



A^/ii(a) + A 



^+^/i2(a) + A^+2/i 



3(a) + ... 



and /ifc(a) are some polynomial conservation law densities. For instance, 
hi (a) = Tja^a^^^~^/2 is a momentum, h2(a) is the Hamiltonian, where the Kodama 
hydrodynamic type system has bi-Hamiltonian structure, and the first of them is 

9h2 



<9^ 



fc = 1,2, 



TV. 



Without lost of generality let us restrict our consideration on the three component case 



w 



Vt 



dJuv) 



Wt 



V 

uw + — 



Thus, the Gibbons equation (see (12) and (36)) is 



op 



Pt- Ox { — 



w 



Then one has 



w{p — u) + 



A 



A„ 



v{p 



u 



Ap, A^ 



{P 



u 



A 



where 



A 



u 



V 



{p — u) —w{p 

Since, the Kodama hydrodynamic type system is homogeneous, then the function A(m, v, w, p) 
must be homogeneous. Then A = {2pdp + 2udu + ?>vd^ + 'iwdyj)\ up to some insufficient 
constant factor (degree of homogeneity). Thus, the equation of the Riemann surface can 
be found in quadratures. For instance, 

Q ^^X = 

^ {p — u) [2p3 — A.up'^ + 2(m^ + w)p 



2uw 



Then (cf. [16]) the equation of the Riemann surface for the Kodama hydrodynamic type 
system is 



X = p 



w 



p — u 2{p — n)2 

Remark: The same procedure can be repeated for any component Kodama hydro- 
dynamic type system. Moreover, suppose the determinant (14) is computed and written 
in the factorized form 



N 



k=l 
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for an arbitrary A^, then the equation of the Riemann surface is given by rational function 
(see [16]) 



A 



11 ''^-^+Z.(p_ai)^- 



(59) 



fc=l ' A:= 

where Bk are some polynomials with respect to flat coordinates a"'. These coefficients Bk 
can be found by substitution the above formula in (11). The corresponding linear system 
is 



N 



rn+l—n 



A; = l,2,...,iV-2, 



n = 2,3,...,N-1, 



dnB'^' = « 

m=n+l 

where B^ = a'^ and 

N~l 

^ a™+^-"a^5^"^ = 0, 57v5" = 0, n = 2,3,...,iV-l. 

m=n+l 

Remark: In the symmetric case (7) all derivatives = (...)Ap can be found imme- 
diately, but in the above case for each A^, one must (step by step) compute all above 
derivatives consequently. If the given hydrodynamic type system is non-symmetric, if the 
generating function of conservation laws (as in the above example) is not given a priori, 
then derivation of the equation of the Riemann surface becomes the very complicated 
computational problem. 

Remark: The equation of the Riemann surface (59) can be written in the totally 
factorized form 

N 



A = (p-ai)i-^J](p-6^'), 



k=i 



where (A^ — l)a^ = S6'^(a). Substituting the Taylor series (17) p' 



(fe) 



Ac'=(b) 



(36) yields the Kodama hydrodynamic type system written in the symmetric form 



.k\2 



2 ' 2{N-1) \^ 



... m 



(60) 



6.2 Cubic Hamiltonian hydrodynamic type system 

The cubic Hamiltonian hydrodynamic type system (24) with the Hamiltonian density 

is equivalent to homogeneous hydrodynamic type system (cf. (60)) 



under the transformation a* = + 2(3T^'^^m^. 
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Following the recipe given above, one can verify that the hydrodynamic type system 
(61) is integrable iff 

""l + 2«E7n ^ ^""l + 2/3E7„ 
and is connected with the equation of the Riemann surface 



A 



2a 



P 



system 



Remark: The mechanical interpretation. Let us consider the classical Hamilton's 

dU dU 



X 



P 



dp dx 

where the Hamiltonian is H = p^ /2 + V {x,t). We seek the Hamiltonian system x = —V^ 
possessing the first integral ("energy constant surface"; see [17]) 



\ = Vq + ViX + V2X^ + ... + Vn-2X^-^ + Vnx'', 



(62) 



where Vk{x,t) are some functions and is arbitrary. Differentiating this equation with 
respect to t yields the hydrodynamic type system (61) (cf. (60)) 



u 



k\2 



.m\2 



U 



2 2(iV+l) 

where the field variables u'^ are coefficients of the polynomial (62) written in the factorized 
form 

A = (x + ^M™) JJ(x-m'=) 
and the "potential energy" is given by the symmetric expression 



V 



E""') +E' 



u 



m\2 



2{N + 1] 

Let us replace (62) on an arbitrary dependence 

A = A(u; x), 

where u^{x,t) are some functions. Differentiating of the above equation with respect to 
one can obtain (see (1), (14); cf. (11)) 

X,^,V{u) = {v1{u)-^^5'^)\u, 

where we use the equation of the Riemann surface A = A(u;/i). In the symmetric case 
(7) the above system reduces to (see (15), (8); cf. (11)) 

\,^.y{M) = A'l{u,^Ji)\k, 

which should he connected with (11) by the substitution (see (15)) 

dip 



dp 
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6.3 The ideal gas dynamics 

The ideal gas dynamics 



= ( Y + ) , vt = d^{uv) (63) 



is connected with the equation of the Riemann surface (see [3]) 



A = + + (64) 



The corresponding Gibbons equation is 

Xt - (/ + u)X^ = ^ 



Pt-d^{ + up 



(65) 



Remark: The equation of the Riemann surface (64) is unique for all values of index 
f3. Indeed, introducing the Riemann invariants 

the ideal gas dynamics (63) can be written in the diagonal form 

1,2 ■'■ r 1 1 2 I / 1 2m 1,2 

rt =-[r +r ± -(r -r )]r^' . 
Then the equation of the Riemann surface (64) in the Riemann invariants is 

A = ^ + ^+^^'-;'^', (66) 



where 



Remark: Two infinite series of conservation laws can be obtained with the aid of 
the Biirmann-Lagrange series (see the next section) at the vicinity of two zeros of the 
equation of the Riemann surface 



_ [z/+(vVl + vV2)2][zy+(vVl-vV2)2] 

4A — , 

u 

while the Kruskal series can be derived at the infinity (A ^ oo, z/ ^ cxd) and at the vicinity 
of another singular point (A oo, u ^ 0). 

6.4 Whitham averaged Sinh-Gordon equation 

A one-phase solution of the Sinh-Gordon 

Uxt = sinhw 
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averaged by the Whitham approach is the two component hydro dynamic type system (see 
[33]) 

r\ = f,\r)rl, (67) 
where the differentials of the quasi-momentum and the quasi-energy, respectively, 

A- < A > 



dp 



v/A(ri-A)(r2-A) 
determine characteristic velocities 

dq{X) 



d\, dq 



i_ < 1 > 



v/A(ri-A)(r2-A) 



dX 



dp{\) 



A=rl.2 



'e(s) 



1 ±1 



where K(s) and E(s) are complete elliptic integrals of the first and second kind, respec- 
tively; = r^/r^ is elliptic module, and 



< a >= — 



adX 



TJ v/A(ri-A)(r2-A)' 



dX 



v/A(ri - A)(r2- A) 



However, this hydrodynamic type system can be enclosed in the framework presented 
in the previous sub-section. Two next theorems can be proved by straightforward calcu- 
lation. 

Theorem 4: The Gibbons equation 



X, 



dq/du dX _ , 
dp/dv op 



connects the averaged (by the Whitham approach) one-phase solution of the Sinh-Gordon 
equation (67) with the equation of the Riemann surface (66), where 



^ 1 2 

p= -+r +r + 



where 



< A > 



In V, 



< A >= 



g = In 



K(£) 



u + (vr^ 



2^2 



< A > 



)2 V' 



In I/, 



Theorem 5: The generating function of commuting flows in the Riemann invariants 



7 /* _|_ 2 P 
z/* - z/(C) 



is connected with the Gibbons equation 

3z/2(A) + 4(ri + r2)z/(A) + (r^ - r^) 



X 



lnz/(C) + 



MA) + 2P)(KA)-KC)) 



A. 



dX 

dp{X) 
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where the potential P of the Egorov metric gu = diP is 



< A >, 



the Egorov metric is 



9ii 



E\s)/K\s) 
2(1 - s2) 
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[1- -E{s)/K{s)f 
2sHl - s2) 



and 



Q(A,C) = 2(ri+r2) + iKA) + ^^^(C) + PlnKA)lnKC) + [r' + r' + ^^(C)]lnKA) 

+ [r^ + r^ + iz/(A)] In u{C) + 2(A + C) ln[i^(A) - z/(C)] - 2Aln u{C) - 2C In i.(A). 



7 Integrable hydro dynamic chains 

The Gibbons equation for the hydrodynamic type system (58) 



op 



p 



2-/3e 



2-pe 



+ ap 



(68) 



where d = a/ {P — Pe + 1) (see (54)), is exactly the same as the Gibbons equation for the 
ideal gas dynamics (65). 
Introducing the moments 

1 



U 



i\(l-l3e){k+l)+l3 



k = 0,1,2. 



{1 - (3e){k + 1) + p 

the hydrodynamic type system (58) can be rewritten as the Kupershmidt hydrodynamic 
chain (see [20]) 

Sf. = + B'B', + [(1 - (3e){k + 1) + (3]B'Bl = 0, 1, 2, ... (69) 

connected with the Gibbons equation (68), where 

Suppose the moments B^ of the Kupershmidt hydrodynamic chain are some functions 
of field variables u^, then TV component hydrodynamic reduction 



(70) 



is an integrable hydrodynamic type system (7) iff the function i?°(u) satisfies some non- 
linear PDE system, which is consequence of the compatibility condition di{dkp) = dk{dip), 
where (cf. (25)) 

-1 



d,B^(n) 



^iy-lSe _ pl-l3e 



/^t -r ^ ^^ky~(3e _ pl-/3e 



(71) 
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If 50(u) = (1 - l3e)T.-i^{u^Y-^^+^/{l3 - /3e + 1), then (71) reduces to (25). The com- 
patibihty condition di{dkp) = dk{dip) creates a nonhnear PDE system on function B^{u) 
only. Its solution is parameterized by arbitrary functions of a single variable. Thus, 
any symmetric hydrodynamic type system (7) can be used for derivation of correspond- 
ing integrable hydrodynamic chain, which have the same Gibbons equation. At the same 
time, all hydrodynamic reductions can be written in a similar symmetric form, but with 
another dependence of r.h.s. functions like B^{u) with respect to field variables u'^. If 
somebody will be able to solve corresponding nonlinear PDE system (which is known in 
the Riemann invariants as the Gibbons-Tsarev system, see [13]), then infinitely many 
symmetric hydrodynamic type systems (7) will be produced. 

Corollary: Let us consider the generating function of conservation laws (see (70)) 



+ —P 



a + i a 



and introduce the parameter a = 1 — (3e. Then + 1 parametric family (A^ parameters 
7;. and (3 for each fixed index a) of hydrodynamic reductions of the hydrodynamic chain 
(69) 

B^, = B^+^ + + [aik + 1) + f3]B''Bl k = 0,l, 2, ... 

is a set of the hydrodynamic type systems (70) 

which are distinct for every value of index /5 (all above hydrodynamic chains are equivalent 
for each fixed index a and for any value of the index j3, see details in [26]). 
Remark: Introducing the moments 

^'-03?^|y^Eo.(«-)'*-"*^ 0.1.2..., 

the hydrodynamic type system (26) can be rewritten as the Kupershmidt hydrodynamic 
chain (see [21]) 

Cf = /3eil + /3C'')-'C^+^ + [i/3e - 1)(A; + 1) +/3]C*^+i[(l +/5C°)-"],, A; = 0, 1, 2, ... (72) 

In this section we proved that hydrodynamic type systems (26) and (58) are com- 
muting flows, then the corresponding [B and C) hydrodynamic chains are commute with 
each other. All other details can be found in [26]. 



8 Hamiltonian chromatography system 

Another generalization of the chromatography system (cf. (23) and [10]) is given by the 
Hamiltonian hydrodynamic type system 

dh 

al = d., — , 2 = 1,2,..., AT, 
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where the Hamiltonian density h(A) and A = TjZk{a^). If this system is diagonahzable 
(see (2)) 

r\ = ^^\v)Tl z = l,2,...,iV, 

then it is integrable (see [6]). Thus, we are looking for corresponding transformation r*(a). 
A direct computation yields (the indexes of the Riemann invariants and characteristic 
velocities /i* are omitted for simplicity below) 



dr 
dzi 



C-Vi 



S ' n 



where Vk{zk) = 472, (Inh')'p(A) = 1/2, <f = p-\A)J:Vkdr/dzk and fi = (W. The Rie- 
mann invariants exist iff the compatibility conditions di{dkr) = dk{dir) are fulfilled, where 
di = d/dzk- Eliminating if and its first derivatives from the compatibility conditions, one 
can obtain the integrability condition 

qj{diqk - dkqi) + qk{djqi - dtq^) + qi{dkqj - d^qu) = 
for every three distinct indexes i,j, k, where qi = {( — ^/)~^- Taking into account 



ViV" V 



the above integrability condition reduces to sole ODE 

VV" = {1 + a)V'' + /3V' + ^ 



(C - v^y 



-1 



(73) 



where Vi = V{zi), zi = z{a}) and p(A) = aA + 5 (a,/?, 7,5 are arbitrary constants). 
Thus, the Hamiltonian chromatography system 



ai = 9,.[h'(A)2'(aO] 



(74) 



is integrable if h = (a = 0), h = In A (a = —1/2) and h = A^ (for all other values a; 
e is an arbitrary constant). 
Then the Gibbons equation 

X, - z"(p)h'(A)A. = Q^iPt- d4z\p)h'{A)]) 

is determined by the equation of the Riemann surface A(a;p), which can be found in 
quadratures 

where ( = z"{p). 

Remark: Introducing the moments 
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the Hamiltonian chromatography system (74) can be rewritten as the Hamiltonian hy- 
dro dynamic chain 

fc+i 

a^ = J2f!:{a)a:, k = 0,1,2,... 

n=0 

determined by the Hamiltonian density h(y4°) (h = , h = InA^, h = (A'^Y) and by 
the Poisson bracket 

where (here we use (73)) 

= [2(1 + a)k + + 2/3fcA^+" + 27^^^=+"-^ 



9 Exceptional (linearly degenerate) case 

In this paper the algebro-geometric approach for integrabihty of hydrodynamic type sys- 
tems is estabhshed. However, this approach is most effective just in case of the symmetric 
hydrodynamic type systems possessing any symmetry operator (see section 6). Then the 
generating function of conservation laws for the symmetric hydrodynamic type systems 
can be found immediately; the integration factor in the computation of the Riemann sur- 
face can be found just if a corresponding hydrodynamic type system is invariant under 
some Lie group of symmetries (like homogeneity). 

Nevertheless, this algebro-geometric approach can be used in all other more general 
and complicated cases, except possibly hydrodynamic type systems which are hydrodynamic 
reductions of linear- degenerate hydrodynamic chains (see [8], [25]). These hydrody- 
namic type systems usually can be written explicitly in the Riemann invariants (see [25]) 
and the conservation law fluxes of their generating functions of conservation laws are 
linear functions with respect to conservation law density p (see [25]) 

Pt = d^[pv{r; X)]. (75) 

For instance, in the limit (3e = 1, the integrable chromatography system (26) 

u\ = eida--, — , , M / -, , i = l,2, ...,A^ 

can be written explicitly in the Riemann invariants (see [10], the formula 5; also [25]; 
the parameters Si and 7^ do not affect on explicit expressions of characteristic velocities 
in the Riemann invariants) 

^* = ^^V^^- i = l,2,...,N. 
Its generating function of conservation laws is (see [25]) 



Pt = ( p- ^ 
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In the same limit (3e = 1, the hydrodynamic type system (58) 



u. 



/e 



in the Riemann invariants is (the parameters 6i and 7^, do not affect on exphcit expressions 
of characteristic velocities in the Riemann invariants) 



r" I rl. 



Its generating function of conservation laws is (see [25]) 



a. 



In this section we consider hydrodynamic type system written explicitly in the Riemann 
invariants 



V r 



1,2,. ..,N, 



such that the characteristic velocities t;*(r) are determined by the unique function v{r;X) 

v\r) = vir;X)\^^,.. (76) 
Then the semi-Hamiltonian criterion (3) reduces to 



diV 



Vi - V 



Vj-V 



(77) 



Suppose this hydrodynamic type system has the generating function of conservation 
laws (75). Then 

A 1 "^'^ 
Oi inp 



Vi - V 



and the semi-Hamiltonian criterion (77) is satisfied automatically. Moreover, the gener- 
ating function of conservation law densities p in such case can be found explicitly 



p = exp 



^ Vk-V 



Thus, any commuting flow 



w r 



has similar generating function of conservation laws 

Pr = d^[pw{r,X)], 

where w*(r) = w{r,\)\x=ri- 

Remark: phase solutions of the nonlinear equations (like KdV, NLS, Sinh-Gordon) 
averaged by the Whitham approach are hydrodynamic type systems belong to the class 
presented in this section (76) 

~ dp 

where the Abelian holomorphic differentials of the quasi-momentum dp and the quasi- 
energy dq are determined on the Riemann surfaces of genus N (see, for instance, details 
in [18]). Thus, corresponding hydrodynamic chains should be linear degenerate. 
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10 Conclusion and outlook 



In this paper we suggest universal (except linearly degenerate case) approach for inte- 
grable (semi-Hamiltonian) symmetric hydrodynamic type systems. This approach is very 
effective because the generating function of conservation laws in such case is unique, this 
is a consequence of the construction ((7) — >■ (9)). In all other cases the first step is a 
computation of the generating function of conservation laws. 

Let us illustrate a complexity of this problem on two component hydrodynamic type 
system (for simplicity and without lost of generality we can restrict our consideration on 
two component case only). Suppose we have nonlinear elasticity equation (this is nothing 
but the ideal gas dynamic written in the Lagrangian coordinates, while (63) is written in 
the Euler coordinates), which is commuting flow to ideal gas dynamics 

Vy = u^, Uy = ^x'-^^Ty- (78) 

This hydrodynamic type system is written in non-symmetric form. Thus, we do not 
know in advance the generating function of conservation laws in this case. Nevertheless, 
obviously, we must seek such generating function in the form 

Py = dxip{u,v,p). (79) 

However, if such generating function will be found, then this generating function is a 
generating function for whole hydrodynamic chain, is not for nonlinear elasticity only. 
Thus, one should seek component hydrodynamic type systems written in the Riemann 
invariants (13) compatible with the above generating function. Thus, we have 

* " dip I dip 

dp lp=P' dp 

Here is a crucial point of the general construction. From the first equation of the 
nonlinear elasticity equation (78) one can obtain (take into account (13)) 

dip 

— \p=pidiV = diU. (80) 



dp 

From the second equation one gets 



dp 



\p=p 



SiU = yf^-^div. (81) 



Thus, we have two choices 



d]P_^lt\ _|_ ^it ^it _|_ ^2— 

ddu dp lp=P* dv o n du dv dp IP=P^ 

iP = dt\_^_dt = — dti — 

dp 'P—P^ dp dp 'P—P^ dp 

These two options are no longer equivalent (in comparison with the symmetric case) 
because we ignore the identity (which is consequence of (80) and (81)) 



dip 
dp 
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which is vahd for ideal gas dynamics only. Since we consider the generating function of 
conservation laws (79) for whole hydrodynamic chain, then depending on this choice, we 
shall be able to construct two different hydrodynamic chains. Thus, for instance, the non- 
linear elasticity equation (78) can be embedded into different Kupershmidt hydrodynamic 
chains (69) and (72). 

Example: The dispersionless limit of the Boussinesq system (see (78), /? = 3; [12], 
[25]) 

Vy = u^, Uy = d^{v'^/2) 

satisfies the Gibbons equation 



2) 



V 



X,. 



dX 
dp 



V 



where ^ = and the equation of the Riemann surface is (cf. (64)) 

,3 



A = /i + 3ti + 



V 



^3) 



Simultaneously, the simplest reduction of the Benney moment chain (see [12]) is again 
the dispersionless limit of the Boussinesq system determined by the Gibbons equation 



where the equation of the Riemann surface is 

A = — /i^ + ?>vji + 3m. 
The substitution ji = —{p + v/p) in this cubic equation yields exactly (83) 

A + 3m + ^. 

Factorizing the cubic polynomial (84) 

A = —p{p — u^){p — u^), 
where p = fl + {u^ + m^)/3, the dispersionless limit of the Boussinesq system (82) 



^4) 



y Q 



satisfies the Gibbons equation 
Ay - (p 



X,. 



dX 
dp 



Py - dx 



-Miu' f - 2uW] 



p^ u-^ + ^ 
~2 3 ^ 



Remark: Generating functions of conservation laws 

,2 



Py = -<9a 



V 



Py = 
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p^ + ^ 



have different sets of conservation law densities H^, which coincide for nonhnear elas- 
ticity equation (82) up to insufficient factors. 

In another paper we present a classification of integrable hydrodynamic chains based 
on the concept of generating functions of conservation laws. For instance, all generating 
functions of conservation laws (79) can be found. Thus, at least two of them are connected 
with the ideal gas dynamics (63); each two component hydrodynamic type system (1) 
must be connected with some function ilj{u,v,p). 

We cannot suggest the recipe how to construct this function ip{u^,u'^, ...,u^;p) for 
any a priori given hydrodynamic type system (1) in general case. However, for any given 
function ip{u^,u'^, ...,u^; p) we are able to reconstruct a corresponding hydrodynamic type 
system (1) together with its commuting fiows. 
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